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Abstract: Our universe may have formed via bubble nucleation in an eternally-inflating 
background. Furthermore, the background may have a compact dimension — the modulus of 
which tunnels out of a metastable minimum during bubble nucleation — which subsequently 
grows to become one of our three large spatial dimensions. Then the reduced symmetry of 
the background is equivalent to anisotropic initial conditions in our bubble universe. We com- 
pute the inflationary spectrum in such a scenario and, as a first step toward understanding 
the effects of anisotropy, project it onto spherical harmonics. The resulting spectrum ex- 
hibits anomalous multipole correlations, their relative amplitude set by the present curvature 
parameter, which appear to extend to arbitrarily large multipole moments. This raises the 
possibility of future detection, if slow-roll inflation does not last too long within our bubble. 
A full understanding of the observational signal must account for the effects of background 
anisotropy on photon free streaming, and is left to future work. 
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1. Introduction 

Inflation is generically eternal [Q, 0. That is, for many scalar field potentials the physical 
volume of inflating spacetime is divergent, with inflation ending only in localized "pockets" 
within the inflating background. This is the case, for instance, when inflation is driven by the 
positive vacuum energy of some metastable "parent" vacuum, in which the vacuum phase of 
our "daughter" universe arises due to tunneling through a potential barrier. The tunneling 
process is described by an instanton that interpolates between parent and daughter vacua, and 
appears as bubble nucleation in the inflating spacetime |^, ^] . Indeed this view of cosmology 
is supported by the present understanding of string theory, which argues for the existence of 
an enormous landscape of such metastable vacua |^, ^, |8|. 

In the standard picture, both the parent and daughter vacua have three large (expanding) 
spatial dimensions. Then the symmetries of the de Sitter parent vacuum suggest the daughter 
bubble should possess a homogeneous and isotropic geometry — in particular it should possess 
an 0(3) rotational symmetry on any homogeneous foliation |^]. Yet the number of large spatial 
dimensions may vary from vacuum to vacuum IC, 11, 12], as is expected in string theory. 
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In particular, one of the large spatial dimensions of the daughter vacuum may be compact in 
the parent vacuum. The size of this dimension can be characterized by a "volume" modulus 
that, during bubble nucleation, tunnels out of a metastable minimum, and subsequently grows 
to very large values. While the resulting bubble is still assumed to reflect the symmetries of 
the parent vacuum, the presence of the compact dimension breaks 0(3) rotational invariance. 
Indeed, we expect the bubble geometry to be toroidal, with 0(2) rotational symmetry in the 
two large spatial dimensions, uniformly wrapping around the compact space. This would 
appear as anisotropic initial conditions to an observer in the daughter vacuum. 

While the formerly-compact spatial dimension remains globally closed in the daughter 
bubble, this will not necessarily be evident to a local observer. In order for the daughter 
vacuum to correspond to our universe, its local evolution should approach an approximately 
0(3) rotationally-symmetric Friedmann-Robinson- Walker (FRW) cosmology, and the circum- 
ference of each large spatial dimension should become and/or remain much larger than the 
Hubble radius. In fact these conditions are easy to satisfy — the latter is akin to solving the 
horizon problem of classical big bang cosmology, and is accomplished by a sufficiently long 
period of slow-roll inflation (after bubble nucleation) along each large spatial dimension [13]. 
(A period of inflation after bubble nucleation is required even in the standard 0(3)-symmetric 
situation, in order to redshift away the large initial spatial curvature of the bubble.) Mean- 
while, it is well known that during such inflation an initially homogeneous but anisotropic 



universe rapidly approaches local isotropy |15|. 



Although the initial anisotropy rapidly redshifts away, background anisotropy present 
at the onset of inflation will generate statistical anisotropy in quantum fluctuations as they 
expand beyond the Hubble radius, and this in turn can modify the spectrum of primordial 



density perturbations |T^, |18|, |19[. During inflation the affected modes are pushed to 
physical scales far beyond the local horizon; however if the duration of inflation is appro- 
priate they will have re-entered in time to form the largest observable scales in the cosmic 
microwave background (CMB). Indeed, seemingly anomalous correlations have already been 



detected among the low-multipole CMB anisotropics pO| , 21, 22, which might indi- 

cate deviations from statistical isotropy in the inflationary spectrum. (It should be noted 
that the significance of these "anomalies" is difficult to assess, and their source(s) could be 
non-cosmological.) In this context a number of models of anisotropy during inflation have 



been proposed; see e.g. [25, 26, 27, 28, M, RQ, 31, 32, However unlike other approaches. 



anisotropic bubble nucleation provides concrete theoretical constraints on the form of initial 
conditions, and serves as a natural extension of the standard inflationary scenario. Further- 
more, it provides an opportunity to conflrm aspects of the landscape/multiverse hypothesis. 

When the background is homogeneous and isotropic, inflationary perturbations generally 
decouple into scalar, vector, and tensor modes (for a review of cosmological perturbation the- 
ory see e.g. [0]). The same cannot be said for fluctuations about anisotropic backgrounds, 
which complicates the corresponding analysis. For this reason we assume that metric pertur- 
bations are suppressed, until the background geometry of the bubble has become essentially 
isotropic, after which standard cosmological perturbation theory can be used. As a concrete 
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model one can imagine the spectrum of isocurvature fluctuations in a sub-dominant scalar 
field, which are much later converted into adiabatic density perturbations (as in the curvaton 
mechanism |35]). Note that this simplification comes at little cost: we are interested not in 
the (model-dependent) amplitude or tilt of the spectrum, but in its statistical anisotropy, and 
this should not depend strongly on the back-reaction of the scalar on the metric. 

Even with this and some other simplifications, the analysis is rather complicated. Al- 
though we obtain an analytic expression for the inflationary power spectrum in terms of an ap- 
propriate set of anisotropic mode functions, we must resort to numerical evaluation to project 
this spectrum onto spherical harmonics. Still, we find that certain patterns are evident: 
whereas the standard (isotropic) picture gives a multipole correlator C^^'mm' = {^e.m^\'m') 
that is diagonal in both H. and ^' and in m and m', and is independent of m, our scenario intro- 
duces off-diagonal components in ^ and I' (when I — ^ = ±2), and introduces m-dependence 
into Cu'mm' (it is still diagonal in m and m'). (These results are not unlike those of [17], 
which studied a Bianchi type I anisotropic cosmology.) Our approximations limit us to a re- 
gion in parameter space where the corrections to Ca^'mm' 

are suppressed relative to the leading 
order terms by roughly the present-day curvature parameter Oj^^j.^. While this greatly con- 
strains the size of these effects, they appear to extend to arbitrarily large giving hope for 
statistically significant future detection |^6|. Note that because the statistically anisotropic 
contributions to Cn'mm' are suppressed by i^curv 

We have focused on statistical anisotropics in inflationary perturbations, however in this 
scenario the background spatial curvature of the bubble is itself anisotropic. In particular, the 
spatial geometry is flat along one direction and open in the two-dimensional planes orthogonal 
to that direction. While our computation of the inflationary spectrum accounts for this spatial 
curvature anisotropy, our projection onto spherical harmonics does not. Indeed, anisotropic 
spatial curvature induces anisotropic expansion, which affects the free streaming of photons 
and thus deforms the surface of last scattering, along with our perception of angular scales 
on it 1 37]. This affects the appearance of the inflationary spectrum, inducing corrections to 
the observed multipole correlator Cwmm'- A full understanding of the observable signatures 
of this model involves combining both of these effects; this is left to future work ]38]. 

The remainder of this paper is organized as follows (an effort has been made to make 
the major sections self-contained). We study the dynamics of anisotropic bubble nucleation, 
within the context of a toy model of modulus stabilization, in Section |2|. The primary goal 
of this section is to obtain the instanton boundary conditions that determine the initial con- 
ditions for the subsequent bubble evolution. However, because our compactification of one 
extra dimension with positive vacuum energy is (to our knowledge) novel, we present the 
model in some detail. In Section ^ we describe the salient features of the post-nucleation, 
background evolution of the bubble, focusing on the (pre-)inflationary geometry (including 
obtaining a simple analytic approximation of the metric). In Section ^ we compute the 
spectrum of inflationary perturbations in a massless scalar field. To better understand the 
observational signatures of this spectrum, we here also perform a basic analysis of its pro- 
jection onto spherical harmonics. Some issues pertaining to the plausibility of observing this 
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scenario are discussed in Section ^. Meanwhile a final summary, including a discussion of 
various avenues for future work, is provided in Section 



Preliminary accounts of this work were presented in |39[] . 
2. Anisotropic bubble nucleation 

We study the possibility that our universe formed via bubble nucleation within some parent 
vacuum in which one of the large spatial dimensions of our vacuum is compact.^ Although 
string theory indicates that our vacuum itself has six or seven compact dimensions, for sim- 
plicity we consider the associated moduli fields to be non-dynamical spectators in all of the 
processes of interest here. Thus the parent vacuum is taken to have one compact dimension, 
the "volume" modulus of which tunnels out of a metastable minimum during bubble nucle- 
ation, which subsequently grows to very large size to create our effectively (3-|-l)-dimensional 
daughter vacuum. Note that the metastable minimum referred to above must have positive 
vacuum energy, so the parent vacuum can decay to our positive-vacuum-energy universe. 

2.1 (2-|-l-|-l)-dimensional modulus stabilization 

We first construct an explicit (toy) model of compactification. The model is based on one 
in 1 10 1 (see also references therein). Our purpose is to demonstrate that there are no basic 



dynamical obstacles to implementing our picture of bubble nucleation, and meanwhile to 
provide a concrete model for future reference. In order to generate a metastable, positive- 
vacuum-energy solution for the volume modulus, we use the winding number of a complex 
scalar field to stabilize the size of the compact dimension. In particular, we consider the 
(3-|-l)-dimensional (hereafter denoted 4d) action 



S = \/~9 



(i? _ 2A) - -K{d^cp*d^^) - - (Iv^p - rj' 



IGttG ' / 2 ' '^'^ 4 



(2.1) 



where g is the determinant of the 4d metric 5^,^, R is the 4d Ricci scalar, A is a cosmological 
constant, and 93 is a complex scalar field for which we allow a non-canonical "kinetic" function 
specified by K. The other terms are constants. Other degrees of freedom, for instance the 
inflaton and the matter fields of the Standard Model, are assumed to be unimportant during 
the tunneling process, and are absorbed into A (and/or g and R). 

The stabilization of the volume modulus of a compact dimension z is most conveniently 
studied using a metric ansatz with line element 

ds^ = e""^ g^h dx^dx^ + L^e* dz^ , (2.2) 



^One might also consider the possibihty that two of the large spatial dimensions of our vacuum are com- 
pactified in the parent vacuum. In this case the volume modulus couples to the Ricci scalar, but because 
the (l+l)-dimensional gravity of the non-compact dimensions is conformally invariant, the theory cannot be 
transformed to the Einstein frame. Although this by itself poses no formal problems to constructing viable 
models, it does complicate the analysis, and so to retain focus we leave this possibility to future work. 
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where ^ represents the modulus field. The effective (2+l)-dimensional (hereafter 3d) metric 
and the modulus ^ are both taken to be independent of z. Meanwhile, the compact 
dimension z is defined using periodic boundary conditions, with < z < 27r, so that it has 
the topology of a circle with physical circumference 27rLe*/^. Note that we have introduced 
the following notation. Any quantity defined explicitly within the effective 3d theory (i.e. the 
theory with the z dimension integrated out), such as the effective 3d metric, is marked with 
an overline. Whereas Greek indices are understood to run over all dimensions, Latin indices 
are understood to run over all but the z dimension. 

We seek a solution for the scalar field ip that stabilizes the modulus ^ with respect to 
small perturbations. The action (|2.1|) gives the equation of motion 



d^[V^K'd>'ip]-V^X^{\ip\^-r]^)=0, (2.3) 

where K' = dK{X)/dX, with X = d^(p*d'^ip. Consider for the moment that is a constant, 
^ = ^p. Then the above equation of motion permits the solution 

where n is an integer, representing the winding number of the phase of if. Ultimately we are 
interested in the dynamics of in which case ^ depends on time and ( [2.4]) is no longer an 
exact solution to the equation of motion. We brush this complication aside by taking 

ri' » ^^'e-* , (2.5) 
so that to leading order we have f = rj e*"^ , which in turn gives 

X = d^ip*d^^ = ^e-* . (2.6) 

Note that we are only interested in dynamics that increase the size of the compact dimension, 
i.e. dynamics that increase ^. Therefore if the above inequality is valid in the parent vacuum, 
it is valid throughout our analysis (we assume that K'(X) contains no poles in X). The 
effective 3d action is then obtained by using this solution to integrate the z dimension out of 
the action. After integrating by parts, this gives 



^ R - - ^ e""'^ - \e-'^'^K{X{'4j)) 



levrG 2 " SvrG 2 



(2.7) 



where we have defined G = G/(27rL), K = 2ttLK, and -0 = ^/a, where a = ViGttG, with 
the rescaled modulus ip being a "canonical" scalar field in the 3d theory. 

Recall that our goal is to find a theory /solution in which the modulus ip is stabilized with 
positive vacuum energy. The satisfaction of these conditions can be verified by studying the 
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effective potential of ip, which corresponds to the last two terms of (2.7). It is sufficient for 
our purposes to identify a single viable model, so let us propose 



K{X) = 27rL {X + K2X^ + k^X^) , (2.8) 

where K2 and K3 are constants. Note that it is inappropriate to consider this form of kinetic 
function as the truncation of a longer series, since we will rely on the various terms being 



comparable in size. In this sense, (|2.8| ) is ad hoc. The choice K{X) = 2/3-^1 + jX — 2/3, 
with /3 and 7 being constants, may seem more natural, since kinetic terms of this form arise 
in certain braneworld scenarios. However we have found that this kinetic term does not 
permit any metastable solutions. Nevertheless, we do not consider this a serious issue. We 
are here focusing on only one very simple model of compactification, as a representative of 
the essential dynamics. It seems reasonable to expect that a thorough investigation of, say, 
the string landscape, would reveal a large number of realistic possibilities. 
The 3d effective potential arising due to our choice of K{X) is 

V{^P) = ^ e--^ + 2.L e"^"^ + ^ e'^"^ + ^ e"^"^) . (2.9) 

Let us now explain the choice of K leading to the above result. If we had used only the 
canonical term, K{X) = X, then the effective potential would have had no stationary points 
at finite ip (recall that we require A > 0). Including the term proportional to K2 creates a 
stationary point for some negative values of K2, but it is always a local maximum. Thus it 
takes two terms in addition to the canonical one to allow for a local minimum. 

The potential V{ip) of ( |2.9D is displayed in Figure for a particular set of parameter 
values: A/Svr = 10"l^ Xq = n^rf/L'^ = 10"i\ K2 = -^o'^ and K3 = O.ieXo"^ (all quantities 
are given in units of G). The value of A is chosen to roughly correspond to the upper- 



bound observational limit of the inflationary potential energy in our universe [40|, while Xq 
is simply chosen to be on the same order as A. The other parameters, namely the coefficients 
that determine K2 and K3, were simply guessed by trial and error. Since n, 77, and L appear 
only in the combination Xq = ri?rf'/Lp', there is no technical problem with satisfying the 
inequality (2.5) for any shape of potential V{il)). Also, we have checked that the above values 



of K2 and K3 satisfy the "hyperbolic condition" of ||4l[ , indicating that the underlying solution 
for ip is stable to small perturbations. 

The effective potential displayed in Figure |l] has a positive vacuum energy local minimum 
at some value ip = ipp. The state with ip = ipp \s the parent vacuum, which appears as 
effectively 3d de Sitter space on scales much larger than 27rLe"^p/^. The daughter vacuum is 
created when ip tunnels through the barrier, to some value ip = ipd, after which tfj accelerates 
from rest and rolls down the potential, with ip ^ 00 as time x'^ — t- 00. 

It is important to this work that the daughter vacuum be created via tunneling (Coleman- 
De Luccia instanton), as opposed to by quantum diffusion over the barrier (Hawking- Moss 
instanton [E2[ ) , because the former gives us a handle on the initial geometry of the daughter 
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Figure 1: The effective potential of the modulus field i}}': see text for details. 



universe. This translates to the requirement that \ V"\ ^ 2i2ttGV at the peak of the potential 
barrier. This condition not hard to satisfy; in particular it is satisfied for the choice of 
parameters given above. Note that this is not the only channel by which the parent vacuum 
may decay — in particular it may transition to any vacuum with winding number n' ^ n 
(see e.g. [10|). Insofar as ip remains trapped in a local minimum, these transitions are of no 
interest, as observers like us could never arise in the resulting 3d daughter vacuum. If such a 
transition eliminates the local minimum or otherwise sets tp free, then we simply consider it 
a more sophisticated version of the case we study here. 



2.2 Tunneling instanton and bubble geometry 

The semi-classical theory of vacuum decay via bubble nucleation is laid out in and we 
find that we can work in direct analogy to that analysis. The tunneling instanton is found by 
studying the Euclidean action, and we begin by constructing a metric ansatz that exploits the 
full symmetry of the parent vacuum. In the 3d effective theory with the compact dimension 
integrated out, this gives the Euclidean line element 

ds^ = df + p\0 [dx' + sm\x) d4>^] , (2.10) 

where ^ is a radial coordinate and x c^nd (j) are angular coordinates on a two-sphere of radius 
p{£,). An effective theory constructed by integrating out a compact dimension is generally valid 
only on scales much larger than that of the compact dimension. In the present case, however, 
the tunneling field is also the modulus of the z dimension. If the 4d manifold factorizes as in 
our ansatz (p.2|), then the modulus ip is independent of z on all scales. Therefore any domain 
wall of must be independent of z, and there can be no topological structures to support any 
z dependence of p in the tunneling instanton. (These statements hold in particular because 
the tunneling field is the modulus of the compact dimension; of course such restrictions would 
not apply to a generic instanton in this background.) 



After integrating out the coordinates x and (and z), the Euchdean action can be written 



S^=A7,j di [-^ {(? + 1) + p2 + 



(2.11) 



where, in the present context, the dot denotes differentiation with respect to ^. The field 
equations of the above action can be written 



^.-8vrG^-^^-yJ (2.12) 

'ii; + 2^ = V'{'^), (2.13) 
P 

where the prime denotes differentiation with respect to if). These are the same as the field 
equations of the standard 4d tunneling instanton, except for some numerical factors Q. In the 
present case, the instanton interpolates between a local excitation of the parent vacuum, in 
which ij) = ipp, and some point on the opposite side of the potential barrier, at which ip = ip^i- 
Inside the daughter vacuum ip continues to grow, but this evolution is not covered by the 
tunneling instanton — the geometry of this region is deduced by matching at the instanton 
boundary, where ip = ipd- As a consequence, the circumference of the compact dimension 
never exceeds IttL e"^'* in the region covered by the instanton. 

It is straightforward to verify the existence of a tunneling instanton by direct numerical 
integration. The boundary conditions are determined by Taylor expanding yo(^) and 

p{0=^-^V{^,)e + --- (2.14) 

m = - ^ ^PdV'i^d) (2.15) 

where the coefficients of the expansions are determined using the equations of motion. The 
value of ijjd is set by trial and error, so that as ^ approaches some value ^maxj the instanton 
smoothly approaches ip{£,) — )■ constant (the excited value just prior to tunneling), and likewise 
pic) ~^ Cmax — Of course the details of the instanton depend on the shape of the modulus 
potential V; to demonstrate the consistency of our model we use the same values of parameters 
as are used to generate Figure ||. The results of this numerical integration are displayed in 
Figure I, and agree with the qualitative description above (note that ^ = corresponds to the 
side of the potential barrier leading into the daughter vacuum). In Appendix^ we analyze 
this tunneling solution from the 4d perspective, substantiating these results. 

The instanton displayed in Figure ^ does not appear to be well-described by the so-called 
thin-wall approximation. Nevertheless, it is not implausible that other modulus potentials 
exist in which the thin-wall approximation is accurate. In such case the analysis can proceed 
in exact analogy to that of Although we do not here propose such a model, we shall assume 
the tunneling instanton can be well-described by the thin-wall approximation to simplify the 
computation of inflationary perturbations in Section]^ (Appendix^. 
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Finally, let us note the geometry of the nucleating bubble . In the 3d effective theory, 
the Lorentzian background geometry is determined by analytic continuation, for instance 
X — + 7r/2 (and C ^ S,^ p ^ which generates the line element 

ds^ = df + /(O [-dx^ + cosh2(x) d(p^] . (2.16) 

This converts the point ^ = into a hypersurface: the future lightcone of the nucleation event 
(at X = 0). The spacetime can be extended inside this lightcone by analytic continuation, 
taking X — >■ X ~ ^7r/2, ^ — )• z^, and p — )■ ip. This creates an open FRW bubble interior. After 
returning the conformal factor e~"^ = e~"^^ of the original metric ansatz ( |2.2[ ), and revealing 
the compact dimension z, the internal bubble geometry has a line element of the form 

ds' = -e-'^mdf + e-"^«)p2(^) ^^^2 ^ sinh2(x) d^^] + L\^^^^^ dz^ . (2.17) 



3. Background bubble evolution 

In the previous section we found that bubble nucleation via modulus decay creates an internal 
bubble geometry corresponding to a (2+l)-dimensional, open FRW manifold crossed with an 
expanding circle, as given by ( 2.17] ). Henceforth we use more traditional geometric variables. 



with respect to which the line element is written 

ds^ = -dT^ + a^(r) [dx^ + sinh2(x) #2] + ^^(r) dz'^ . (3.1) 



We now see that the metric is of the Kantowski-Sachs form (case 2 of ^^), a Bianchi type 
III homogeneous but anisotropic geometry. 

To understand the internal dynamics of bubbles like ours, we add to the Lagrangian some 
matter degrees of freedom. These, we presume, have stress-energy tensor components that 
are all initially negligible next to those of the inflaton. Furthermore, we presume the inflaton 
field (initially) reflects the symmetries of the tunneling instanton; in particular we take the 
inflaton to be static at the instanton boundary and independent of the initially-compact 
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dimension z. It then suffices to treat the matter as a perfect fluid, with stress-energy tensor 
T^v = diag{— /9, p, p}, where the energy density p and pressure density p are (neglecting 
for the moment quantum fluctuations) functions of time alone. The field equations are then 



a ah 1 



a b 



Vib) + p 



b d ( a b 
T + ■ 



a \ a 



P+(2^ + ^)(P + P) 



SnGbV'ib) 
0, 



(3.2) 
(3.3) 

(3.4) 



where here and below the prime denotes differentiation with respect to b, the dot differen- 
tiation with respect to bubble proper time r, while the term V{b) is described below. In 



the second equation, (3^), we have combined two of the Einstein field equations in order to 
express a relation that will be useful later. 

The anisotropy of the bubble is encoded in the initial conditions for bubble evolution, 
which are in turn established by the boundary conditions of the tunneling instanton. These 
set the spatial and time derivatives of all non-metric fields to zero, so at r = we have 



Pinf 



p = 0, 



P 



(3.5) 



where as noted above we assume the initial matter density is dominated by the inflaton. Be- 
cause the "volume" modulus of the z dimension is the tunneling field, the instanton boundary 
conditions require that 6(r) approach a constant as r — >• 0, i.e. 



0, 



(3.6) 



Meanwhile, referring to ( 2.14| ) and ( p. 15 ) and being careful to track the various variable 
redefinitions, we find that as r — )■ 0, 



0. 



1. 



(3.7) 



In Section |^ it was convenient to absorb all of the matter degrees of freedom into the 
cosmological constant A. Moreover, in the 3d effective theory with the z coordinate integrated 
out, A coupled with the modulus ip in the 3d effective potential V{tlj). Now we have returned 
to the 4d picture, in which A is a bare cosmological constant, and have introduced p and p 
to generically characterize the matter content of the universe. Therefore it is convenient to 
shift the cosmological constant A into those quantities. The other ingredients of the (toy) 



compactification model of Section 2.1 generate a 4d "potential" for b 



(3.8) 



where we note the change of variable dr = e "^/^ dx^, which also affects the form of V{b). 
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For simplicity we treat the inflaton energy density as equivalent to vacuum energy density, 
Pinf = constant. (Thus, according to the above assumptions, pi^i is almost identical to A, 
the difference being the present dark energy.) It is convenient to also assume the inflaton 
energy density dominates over the modulus potential, even right after bubble nucleation, i.e. 
Anf S> V{bd). This allows us to obtain the simple analytic solution, 

"('^) = ^hil sinh(Finf r) , 6(t) = ba cosh(Finf r) , (3.9) 

where H?^^ = 87rG/9inf/3. (Note that the inequality pinf ^ V{bd) is satisfied by the choice of 
parameters used to generate Figure ||.) We have chosen the normalization of the scale factor 
6(r) so as to conform to the initial condition given in (^.6|). 

One might notice that this solution resembles a slicing of 4d de Sitter space. In fact, in 
the (4+l)-dimensional Minkowski (ds"^ = —dt^ + du^ + dv^ + dx^ + dy^) embedding: 

t = HrJ sinh(/7inf r) cosh(x) (3.10) 

= -f^inf COsh(i?inf r) COs(6dFinf z) (3.11) 

^ = -^inf cosh(i?inf r) sm{bdHin{ z) (3.12) 
X = i^j^£ sinh(i7inf r) sinh(x) cos((?i)) (3.13) 
y = ^hd sinh(i?inf r) sinh(x) sin((/)) , (3.14) 

the bubble coordinates {r, 4>^ sit on the hyperboloid —t"^ + + v"^ + x"^ + y"^ = H^^, on 
which (|3.9D gives the scale factors of the induced metric. However this is merely a curiosity; 
it arises only because we take the limit ^nf ^» ^(^d), with pinf = const. Indeed, closer 
inspection reveals that that the Minkowski coordinates u and v are not periodic with the 
bubble coordinate z (due to the factors of ^^-ffinf), indicating that the bubble geometry covers 
only a subset of the hyperboloid, with periodic boundary conditions.^ 

Note that the bubble geometry features two forms of anisotropy: expansion anisotropy — 
stemming from the two distinct scale factors a and b — and spatial curvature anisotropy — 
stemming from the {x, </>) plane being open while the orthogonal z direction is flat. During 
inflation both forms of anisotropy rapidly redshift away. Although the presence of anisotropy 
makes choosing the definition of a scalar curvature parameter somewhat ambiguous, we find 
it convenient to write ^ 

^curv = ^g2jj2 ' (3.16) 

where H = a/a and, following the assumptions above, here and below treat the effect of the 
modulus potential as negligible next to the above term. Notice that the relationship between 



^We can also define the time coordinate ^ = H-^^ sinh(_H'inf r), in which case the line element becomes 

ds^ = - (1 + Hltey' de + e [^X' + s.mh\x) d<t>^] +bl{l + Hltf) dz" , (3.15) 

which is also seen as 4d de Sitter space when = -ff~f . Notice that we can now read ofT the 3d dynamical 
fields of ( IITtI ): p(e) = e (1 + Hlt^)^^'^ and = ""^ ln(l + Hlf^), with L = ba- 
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^curv and l/a^H^ involves a factor of 1/3 that is not present in the isotropic case. Meanwhile, 
the expansion anisotropy can be parametrized by 

ah , ^ 

h^--T. 3.17 

a 

In (^^) we have combined two of the Einstein field equations so that we can now write 



h+[2^ + ^^h = \. (3.18) 
o b I 



Thus we see that anisotropic spatial curvature sources anisotropic expansion. 

While it would be interesting to understand any observable effects of late-time spatial 
curvature and expansion anisotropy (some of which have been explored in 47, 37|), 



for the purpose of this paper we set such questions aside. Indeed, to gain some understanding 
of the statistical anisotropy of inflationary perturbations in Section we project the spectrum 
onto a two-sphere "of last-scattering," treating the background evolution as if the expansion 



anisotropy h can be ignored. As suggested above and shown in greater detail in [^, this 
approach is not entirely self-consistent, which is to say the actual surface of last scattering is 
deformed by the presence of anisotropic expansion h. However the analysis allows us to gain 
a qualitative understanding of some of the effects of anisotropic bubble nucleation; we leave 
the more complicated, complete treatment to future work. 

For the moment, we simply note two important relationships for later reference. The first 
is the leading-order comoving distance to the surface of last scattering. 



/ -^^^—Tr^(^-^V^<0.50, (3.19) 
Jr, a{T) aoHo 



where, to obtain the leading-order expression, we have ignored the effects of spatial curva- 
ture and background anisotropy on the evolution of a (we have included the effect of dark 
energy, modeled as cosmological constant with i^A/Qm = 2.85). Here a subscript "0" desig- 
nates quantities evaluated at the point of present detection while a subscript designates 
quantities evaluated at the surface of last scattering, and the last inequality follows from the 
observational limit ^1°^,^ < 6.6 x 10"^ (WMAP+BAO+SN with w = -1 prior, treating the 
95% confidence level as if it were a hard bound) however it should be emphasized that 
late-time anisotropic expansion induces a quadrupole in the CMB, itself appearing at order 
which absent cancellations provides a stronger constraint on ^^utv 137). 



The second important relationship comes from realizing that periodicity in one of our 
large spatial dimensions has not been observed (see e.g. |Q), indicating that the physical 
circumference of the closed z dimension must be greater than the physical diameter of the 
surface of last scattering, or 

dr 



27TK>2a^ —-. (3.20) 
a(r) 
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Ignoring the subleading effect of late-time anisotropic expansion, we see from ( p.9| ) that we 
can write 6^ = bdH^^i o-*- Inserting the result of the last paragraph then gives the constraint 

O[!,„<0.27(6di/inf)'. (3.21) 

Thus, as the quantity ft^-f^inf is decreased, the maximum allowed present-day curvature pa- 
rameter is decreased, along with the size of all related effects. However it should be noted 
that this result depends crucially on the "redshift factor" hdH\^f between the two scale factors 
a and 6, which in turn relies on our limit of approximation to obtain the analytic solution 
( [3.71 ). It is not clear how this constraint is modified outside of the limit /Ojnf ^» V{bd)- 



4. Inflationary perturbations 

4.1 Power spectrum of a massless scalar 

We would like to avoid the formidable task of developing cosmological perturbation theory 
about the anisotropic background of ( p.l[ ). Therefore, we assume that metric perturbations 
are negligible, until a period when the local background anisotropy is negligible, during which 
standard cosmological perturbation theory can be used. This could very well be a generic 
feature of realistic models of inflation, since at least in isotropic models the metric perturba- 
tions are suppressed (by the flrst slow-roll parameter) relative to inflaton fluctuations, until 



their wavelengths have expanded far beyond the inflationary Hubble radius (see e.g. |48|). 
However to be concrete and self-consistent we focus on the fluctuations in a subdominant 
scalar fleld a. The fluctuations in a may be converted into adiabatic density perturbations 
by a variety of proposed mechanisms; for instance a could be a curvaton Because we 
ignore metric perturbations, we must set aside the interesting question of what modiflcations 
to the spectrum of tensor perturbations one might observe in this scenario. 

For simplicity we assume that during inflation within the bubble the mass of a is negligible 
[m^ <^ Hi^f), as are the interactions between a and any other matter flelds, so that the only 
terms involving a that appear in the Lagrangian are given by 

S„ = - j^<fx ^g^-'d^ad.a . (4.1) 



The background metric Qf^i, is given by (3J.); however we flnd it convenient here to work in 
terms of the "conformal" time r]c, where 

Note that we use the term "conformal" loosely, as this time parametrization does not make 
the metric entirely conformally flat. Instead, we have the line element 

ds'^ = a^{r^^) [-di]l + dx^ + sinh2(x) d<j?] + }?{r]c) dz^ . (4.3) 
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We solve the equation of motion of a by separation of variables, writing 



ab y ah 

where the placement of indices anticipates results below. The factor (a6)~^/^ is introduced so 
that the temporal mode functions T^^ are canonically normalized. The resulting equations 
for $r and Zg are very simple: 

^';{<P) = -r2 $,(^) , Z'^{z) = -s" Zs{z) , (4.5) 

where here and below we use a prime to denote differentiation with respect to the lone 
argument of the function that it accents. The factors r"^ and parametrize the separation 
of variables. A complete set of orthonormal solutions are 

M^) = , U^) = , (4.6) 



where the periodic boundary conditions on (j) and z dictate that r and s must be (positive or 
negative) integers (recall that we scale the z direction so that it is periodic on < z < 2tt). 
Meanwhile, the differential equation for the Xq^ can be written 

(1 - c2) X;,(c) - 2cX'q,{c) = + + Xqric) , (4.7) 

where the sum 1/4: + corresponds to the third separation constant in the above separation 
of variables. We have defined the variable c = cosh(x) so as to obtain the Legendre equation 
(as per our convention the primes now denote derivatives with respect to c), the solutions of 



which are the associated Legendre functions P[g_i/2 and <5[g_i/2 tMI- interest in the 



solutions P[g_i/2^ with q real and positive, which are finite and stationary in the limit c — )• 1. 
A convenient normalization is [5C] 

XM = ^^^^^P[q_,/2{c), (4.8) 
where F denotes the gamma function. The Xg^ then satisfy the orthonormality condition, 

dcXgr{c)X;,^{c)=S{q-q'), (4.9) 



poo 



where 5 denotes the Dirac delta function. Functions of this sort are studied in [50|, and it 
is straightforward to extend those analyses to the present situation. Using the shorthand 
C/gr-s(x) to denote the combined spatial mode functions, one can show that 

Uqrsi^,) f/;.(x,) = ^^^^ Y: , (4.10) 
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where H12 = cosh(xi) cosh(x2) — sinh(xi) sinh(x2) cos((/)i — 4>2)- By following the analogous 
calculation in [^] , this "addition theorem" can be used to prove the completeness relation 

/ dqS^ Uqrs{ci, (t>l, Zi) U*j.g{c2, h, ^2) = S{ci - C2) (5((/)l - 4>2) S{zi - Z2) ■ (4.11) 

The time evolution of a is given by the mode functions T^^, which satisfy 

= 



oo? 1 2 1 ^ lab lb la? 1 &^ i / \ 



(4.12) 



To proceed we must specify the (conformal) time dependence of the scale factors a and 
b. For simplicity we adopt the analytic solution of Section ^, given by ( p. 9] ), which treats 
the inflaton energy density as constant. In terms of this solution, the conformal time is 
r/c = In [tanh(ffinf r/2)], which runs from minus infinity to zero as the bubble proper time r 
runs from zero to infinity. The two scale factors are then given by 



«(^c) = --f^inf csch(7?c) , 6(r/c) = -bd coth(?7c) • 



(4.13) 



Some simplification occurs when we plug this solution into (4.12), and in the end we find 

1 



T"(r/c 



- 2 csch^ (r?c) + ( /^^ - ^ ) sech^(r/c) 



t:;.(??, 



qsVIc) , 



(4.14) 



where we have defined n = s/bdH[^f. The solutions to this equation can be written 



4-^-2' -4-y+2'^-^^' '"'^ 



-2*''C2" coth(?7c) sech*«(r/c) F 
C^iTi{r,,) + Cl^T2{r]c), 



'4 + y-2'-4 + Y+2'^ + '^' '"'^ ^"^^ 



(4.15) 
(4.16) 



where F = 2F1 denotes the hypergeometric function I^Oj. The index v is explained below. 
We have introduced the shorthand notations J-i and 7^2 for later convenience, and have chosen 
to represent the solutions so that, in the limit ?7c — >■ — 00, 



For later reference we also note the asymptotic behavior as rjc — s- 0, 



a{r]c)b{r]c) 
^ a{r]c)b{'qc) 



y'?r(|)r(i-ig) 
hd r(f-f + f)r(f-f -f) 
R-, 2--gr(|)r(i + ig) 
bd r(f + f + ^)r(f + f -f) 



(4.17) 

(4.18) 
(4.19) 
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The initial conditions for the inflationary perturbations are set by the choice of integration 
constants C|' and • Different choices of initial conditions correspond to different choices 
of the quantum vacuum state l5l| — a particularly attractive choice is that of the Bunch- 



Davies vacuum |52], which maps de Sitter mode functions onto zero-occupation Minkowski 
wavefunctions in the small-scale limit. Within the context of bubble nucleation in an inflating 
background, the Bunch-Davies state is determined by tracing the evolution of mode functions 
back into the parent vacuum and Klein-Gordon normalizing the positive- frequency modes on 
a Cauchy surface E^. The analysis is rather tedious, and we relegate it to Appendix 



y. The results are inserted at an appropriate point below. 

Having solved for the mode functions of a, to obtain the spectrum of inflationary per- 
turbations we promote o" to a quantum operator, a, and compute the two-point correlation 
function. The analysis can proceed in exact analogy to the standard formalism (see e.g. p^ ). 
In particular, we can express (T as a mode expansion of creation and annihilation operators, 

1 



o-{Vc, x) 



(4.20) 



where a^qrs and dtqrs satisfy the appropriate analogues of the standard commutation relations. 



[(Ivqrs, (''v'q'r's'] 



d{q- q ) 6rr' 5ss' 



i^vqrs^ ^v'q'r's'1 



(4.21) 



as do the canonical field operator yaho and its conjugate momentum field operator vr. The 
index v is explained in Appendix]^; it takes one of two values, which can be denoted ±. Prom 
( [4.201 ) we see that the "Fourier" transform of <5"(?7c, -2) can be written 

(T(r/c, 9, s) = , ^ ^ = \^\sknc) avqrs 



iVc) al 



v{-q){-r){-s) 



a{'nc)b{ric) 

which gives the equal-time momentum-space two-point correlation function 



(4.22) 



{a{r]c, q, r, s)a''{r]c, q' , r', s')) 



-6{q- q')5rr' Sg 



(4.23) 



a{ric)b{'nc 

Each "Fourier" mode Uqrs associates the set of separation constants {g, r, s} with a set of 
comoving distance scales in the X; 4', and z directions, according to the characteristic scales of 
variation of Uqrs{x^ 4'^ After these comoving scales grow larger than the Hubble radius, the 
mode "amplitude" Ylv l^gsP/*^^ rapidly asymptotes to an r/c-independent constant. Because 
of this, and because all observable scales in our universe first expanded beyond the Hubble 
radius deep in the inflationary epoch, for practical purposes we can safely evaluate the two- 
point correlator ( 4.23| ) in the limit rjc — )• 0. This gives the power spectrum 

^ ^\rU^c)? (4.24) 



Pr 



qs 



r)l™o a{rjc)b{7]c) 
8 bd sinh^(7rg) 



vr cosh(7rg) 



r( + ^ 

4 2 2 



^ / 5 iq 
r - + — 
4 2 
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COS (vr \/l — Re 



22»g r(i -iq- yr^) r(i -iq+ v^w^) 
1^2(5 _ il _ l±\ r2/'5 _ H _i_ 

^ \A 2 2/ V4 2~'~2/ 



(4.25) 



H? 



2 bdHmf 



where we have inserted the results from Appendix |B|. The last relation corresponds to taking 
the asymptotic limit of large q and fx.^ 

To begin to understand this result, first note that the separation constant q characterizes 
comoving radial distance scales in the {x, (p) plane, and thus relates to the standard Cartesian 
wavenumbers and ky via (7^ ~ A;^ + ky in the flat space (large q) limit. Meanwhile, physical 
scales in the z direction are redshifted relative to those in the {x, <p) plane due to the dissimilar 
evolution of the scale factors a and b at early times. The late-time effect is an additional 
factor of b{r]c) / a{r]c)\rjc-^o = bdHi^f relating physical and comoving distances; therefore we 
relate the separation constant s to the late-time Cartesian comoving wavenumber kz via 
s ~ bdHi^ikz- Thus the term in parentheses in ( 4.261) approaches (fc^ + ky + A;^)~^/^, which 
gives the standard scale-invariant power spectrum (the factor of \/bdH\^f will cancel when 
integrating over wavenumbers to compute observables in the isotropic limit). 

Although we have just remarked on the congruence between the asymptotic limit of our 
result, ( |4.26| ), and the result from standard, flat, isotropic inflation, the two inflationary spec- 



tra are not the same. The power spectrum ( 4.26 ) is expressed in terms of the anisotropic 
"Fourier" modes C/grs(x), while for instance the isotropic Cartesian Fourier modes are e**^'^, 
where y = {yi, y2-, ys} are Cartesian coordinates. This difference affects the observed spec- 
trum because the (lack of) correlations implied by the delta functions 8{q — q') Srr' 5ss' in the 
two-point correlator ( [4.23| ) are different than the (lack of) correlations implied by (5*^'^^(k — k') 
in the standard isotropic case. 

4.2 Projection onto spherical harmonics 

At present the spectrum of primordial density perturbations is most tightly constrained by 
measurements of fluctuations in the temperature of photons streaming from the surface of last 
scattering, i.e. CMB fluctuations. It is conventional to use spherical harmonics to cover a two- 
sphere representing our fleld of vision (here parametrized by angular coordinates < 9 < it 



^Instead of performing the involved analysis of Appendix |b|, one might determine the integration constants 
Ci and C2 by studying the behavior of T^^ at very early bubble times, rjc — >■ —00, and on very small bubble 
scales, q 00, and equating it with that of a free field in Minkowski space, T^^ — >■ {2q)~^^^e~^''^'' . (Note that 
the equation specifying T^g takes the form T^^ = — T^^ at early conformal times — hence the "Minkowksi" 
wavefunction is independent of s.) This corresponds to choosing the analogue of the so-called "conformal" 
vacuum, see e.g. It has the significant drawback of predicting an energy density of fluctuations that 

diverges as r^c — > —00, thus converting the mere coordinate singularity at the instanton boundary into a physical 
singularity. Proceeding nevertheless, in this situation one finds no cause for the index v, and inspecting the 



asymptotic behavior of and T2 in (4.17) it is apparent that one would choose Ci = (217)-^/^ and C2 = 0. 



The resulting power spectrum is given by the first term of (4.25), but divided by coth(7rg). It approaches the 



same limiting behavior as the Bunch-Davies vacuum, (4.26) 
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and < (/) < 2tt). The observables are therefore the multipole moments 

aim = JdCd<p y4(C, </)) 5t(C, <t>) , (4.27) 

where St represents the temperature fluctuation and for convenience we have defined C = 
cos(0). The (orthonormal) spherical harmonics are given by 



The temperature fluctuation 6t is sourced by density perturbations on the surface of last 
scattering, but the observed CMB perturbations include secondary effects, which are incurred 
as photons stream from last scattering to the point of detection. For simplicity we here ignore 
these effects. In fact, we completely ignore the presence of late-time anisotropic expansion 
h = d/a — b/b. In this case, null geodesies radiating from the origin of coordinates see a flat 
metric, equivalent to that of a cylindrical coordinate system. Meanwhile the corresponding 
flat cylindrical coordinates can be related to flat spherical coordinates, in terms of which the 
surface of last scattering has fixed radius q-^ (again, neglecting h). Thus we obtain 

X.iC, cP) = gWl-e, MC,<P)=^, = ' (4-29) 

which reflects a particular choice of matching between the cylindrical and the spherical coor- 
dinates (recall that C = cos{9)). The factor l/bdH[^f in z comes from the ratio of scale factors 
a and b when matching onto the isotropic cylindrical coordinates. 

We emphasize that while this procedure allows to create a picture of the inflationary spec- 
trum in terms of spherical harmonics, the actual observed spectrum will contain corrections, 
coming from the anisotropic expansion between the surface of last scattering and the point 
of present detection. In particular, the presence of anisotropic expansion deforms the surface 
of last scattering away from the surface defined by ( [4.29 ), and perturbs the trajectories of 



geodesies as they radiate away from the point of observation. Our analysis can be viewed as 
a first step toward understanding the inflationary spectrum; a more complete analysis of the 
observational signatures is left to future work. 

To be precise, we simply take the temperature fluctuations to be given by 

SriC, </>) cc a{iic = 0, x.(C, 4>)) , (4.30) 



where = {q^,, (j)} designates the coordinates of the surface defined by ( 4.29 ). Here we 
have assumed that the isocurvature fiuctuations in the light scalar field a ultimately directly 
source the adiabatic CMB perturbations, and have ignored a model-dependent proportionality 
constant. Also, in addition to ignoring the expansion anisotropy h, we have ignored the 
evolution of primordial perturbations after they enter the Hubble radius. 
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We are actually interested in the multipole correlation function Cuimm' = (^^m^l'm')- 



Expanding a{r]c. = 0, x^(C, (p)) in terms of mode functions Uqrs as in ( 4.20| ), we obtain 

C«'mm' = / dCi d(2 #1 #2 dq ^ Pqs >£m(Cl, h) ^<?'m'(C2, 4>2) Uqrs{Cl, 4>l) f^grs(C2 



where Pqs is the power spectrum, given by ( pHD , and 

r(i + ig' - r) 



27rT{iq) 



P 



iq-l/2 



coshy^ q'^ (1 — C^) 



Jr<f>+ifi 



(4.31) 
(4.32) 



where fi = s/bdH\ai- Plugging the and C/g^s into ( 4.31 ), we can immediately perform the 
integrations over and 02- This gives 



mm' — '-'mm' 



{2e + l)(2f + !){£- m)!(f - my. 
167r2 (^ + m)!(f + m)! 



r(i +iq - m) 



pm 



T{iq) 

xPr*(Ci)^™(C2)e^^^*(^^~^2). 



COshJ^2(l_^2) 



pm* 



COs\ijQl{l-Cl] 



(4.33) 



Because of the complexity of this result, we here limit our attention to two basic goals: 
demonstrating that, as expected, the Cwmm' approach the standard, isotropic results in the 
limit of many e-folds of inflation, and understanding the qualitative features of the low- 
multipole Cwmm' via approximate numerical integrations. 

The physical distance to the surface of last scattering is fixed by the late-time big bang 
expansion history. To model this, we treat the physical distance from the origin to the surface 
defined by ( 4.29| ) as fixed, a*^)* = constant, so that g^, decreases as the duration of slow-roll 
inflation increases. Thus, the limit of long-duration inflation corresponds to the limit of very 
small Qi,. In this limit only modes with large wavenumbers q are relevant to observation, and 
the normalized Legendre functions approach normalized Bessel functions, 



^(| + ^g-r) 
T{iq) 



pr 

^iq-1/2 



coshy^ £.2 (1 _ ^2) ^ y/qJrig^qV'^ - C) , 



(4.34) 



where to be clear Jr is the rth-order Bessel function of the first kind |49], and we ignore the 
unimportant overall phase. (This result is most easily obtained by taking the appropriate 
limit of the underlying differential equation.) With this substitution the angular integrations 



over and C2 can be performed analytically, using a convenient mathematical equality [56|: 



dCe'y-°<-)<P[^{C)J^ 7sin(a)v^W 



2i 



l—m Tim 



Pr[cos(a)]i,(7), 



(4.35) 



where denotes the spherical Bessel function of order I This result holds for positive 
and negative m, and any < a < tt. Since P™* = (— l)"™'^™, this equality also holds if we 
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replace the Legendre polynomial with its complex conjugate on both sides. Note that it is 



always possible to choose the a and 7 of ( 4.35 ) so that 7cos(a) = ±iJ,Qi, and 7sin(a) = qg^ 
in particular one sets 



7 = Q^\f(f^r[P' , cos(a) = , ^ . (4.36) 



Putting all of this together, we obtain 



\2l + l)(2f + 1)(^ - m)!(f - m)\ 
7r2 (^ + TO)!(f + m 



mm' ^mm' 



s 



X i. Fr (^^) PT (^^) . (4.37) 

Our next approximation is to replace the sum over integers s with an integral over real 
fi = s/bdHi^f. Demonstrating the strict validity of this approximation is tedious; however 
intuitively we expect it to be accurate at least insofar as q-^ <^ bdH^ai- This is because modes 
with wavelength A ~ 1/s should not contribute significantly to observables on scales g <^ \, 
and at the same time the discrete spectrum should be well-approximated by a continuum 
when s » 1, or /X » l/hdH-^^^. Converting the sum over s into an integral allows us to make 
use of a convenient variable transformation, defining k and Q according to 

q = (k/g^) sin(G) , ^ = (k/g^) cos(G) , (4.38) 

where k > and < < vr. The Jacobian of the transformation gives dq dfi = g^'^ k dk dQ = 
g~'^ [k/ sin(0)] dkdcos{Q), so that Cee'mm' can be written 



^ ^ /(2£ + l)(2f + l)(^-m)!(f-m)! ^ /■ „ , 

C«'mm' = 6mm' \l- ^2 ^ ^) , ^ ^) I ^<^^-f J dk dcOs{Q) Pqs{k, G) 



X h{k)k'{k) Pr[cos(G)] P7[cos(G)] . (4.39) 

After this variable redefinition one can also show that in the limit of small g^,, the power 
spectrum approaches Pqs{k, G) — >• {Hin{/2bd) gl/k^ . The integrals over cos(G) and k can 
then be performed, giving the standard flat, isotropic result: 

Now let us turn to computing soniG of the low-multipole components of Cn'mm' • 

The 

direct numerical estimation of ( [4.33| ) converges very slowly, given our limited computational 
resources. Nevertheless, we can proceed as above, replacing the sum over s with an integral 
over ^. This by itself does not improve the convergence of the numerical evaluation, however 
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£ = 2 


e = 's 


£ = A 


£^5 


£ = 6 


£^7 


^ = 8 


£ = 9 


I' = £ 


-1.2 


-1.1 


-1.0 


-1.0 


-1.0 


-1.0 


-1.0 


-1.0 




























£' ^e±2 


0.33 


0.35 


0.36 


0.37 


0.38 


0.40 


0.40 


0.40 



Table 1: The multipole correlator contrast SCu'mm' for several values of £ and £' , ni = m' — 0, in 
units of ilcurv, using — 0.01 (see the main text for details). 



if we now replace the power spectrum Pgs of ( 4.25| ) with its asymptotic limit ( 4.26| ), the 
integral over /i can be performed analytically, 



2q.\Ci - C2I 



Ki ig.qlCi - C2I 



(4.41) 



where Ki is the (first order) modified Bessel function of the second kind. 

As explained above, we expect replacing the sum over s with an integral over fi to be 
accurate when <C b^Hi^f. Meanwhile these quantities are observationally constrained by 
( pj| ) and (|^, which comb ine to give < 3.2 bdH^^f. Thus there is an interesting region 
in the parameter space that cannot be probed by this approximation. On the other hand, 
the asymptotic limit of the power spectrum Pgs, ( 4.26| ), is approached rapidly as q increases; 
see Figure ^. We have computed a few elements of Cwmm' without this approximation, using 
bdHini = 1 and = 0.5, and the two agree at about the 10% level. Note that, for fixed 
bdHi^f, the accuracy of these approximations improves as g^, is decreased. 

Putting everything together, our final expression for the Cu'nim' is 



>t2t + l)(2f + ™)!(f - m)! f ^ 



167r2 {£ + m)\{l' + m)\ 



r(^ + iq — m) 



prn 

^iq-1/2 



T{iq) 

xPr{Ci)PF{C2)Ki{e,q\Ci-C2 



COsh./^,2(l-C?) 



pm* 



COshj£,2(l_^2. 



(4.42) 



To understand the deviations from isotropy, it is convenient to define the quantity 



(0) 

ii'mm' 



which gives the correction to the "isotropic background" c]^^^^ of ( [4.40| ) , in units of a relevant 
component of C^^jnm- Note that within the context of the above approximations, the only 
observable that enters dCwmm' is the comoving distance g-^. Comparing numerical results 
for several different values of we find (roughly) SC^i'mm' gt when g^, < I. Meanwhile 



(4.43) 



( |3.ig ) gives gl ~ STri^urvJ thus we conclude the non-zero components of dCwmm.' scale with 
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Figure 3: Non-zero values of the the multipole eorrelator contrast SCw; 
units of ilcurv, using g^, = 0.01 (see the main text for details). 



for 



3, 6, in 



the present day curvature parameter r^curv Table |l] we display SCu'mm'^ in units of J^iurv, 
for several values of I and choosing for simplicity m' = m = 0, and setting = 0.01. 

The most salient feature of Table || is the existence of off-diagonal terms with respect 
to i and i' . This is qualitatively not unlike the results of [17|, in which the phenomenology 
of a Bianchi type I anisotropic universe was studied. As was the case in that scenario, we 
find that the off-diagonal terms of Cee'mm' do not fall off very rapidly relative to relevant 
components of the background cj!^}^^,. Indeed, the sequence of dCu'mm' with £' = i ± 2 
appears to increase with increasing £. Although computing 5Cwmm' at large ^' = ^ + 2 is 
extremely time consuming, we have explored several values of £' up to i' = 20 and observed 
that the sequence asymptotes toward a small constant value, about 0.4 (in units of riSrv)- 
We have not displayed components of SCwmm' with £' — £ > 1 and odd, which we expect, like 
£' — £ = 1, to give precisely zero. The components with £' — £ > 2 and even feature very slow 
numerical convergence, but also appear consistent with zero. 

The failure of off-diagonal terms of SCei'mm' to approach zero in the limit of large £ may 
at first seem to conflict with the usual interpretation that increasing multipole moments £ 
correspond to probing smaller physical scales. To understand this effect, we must recall that 
we are evaluating the anisotropic modes on a surface that sits at fixed, non-zero q-^. The 
asymptotic behavior of ( |4.34D , in which the anisotropic modes converge to isotropic ones, 
corresponds to the limit q ^ oo, keeping g^^q constant, which must not be confused with 
the limit q — t- oo, keeping g^ constant. In the latter case, there is always a region on the 
two-sphere (corresponding to large intersection with the open plane) over which there is a 
discrepancy between the anisotropic and isotropic modes, even in the limit q — )• oo. Since 
the spherical harmonics receive support over the entire two-sphere (regardless of £), this 
discrepancy corresponds to statistical anisotropy in the projection onto spherical harmonics. 
We emphasize that the above comments concern the ratio that appears in dCu'mm'] the 
observed multipole moments, Cwmm'^ indeed decrease at increasing £ and £'. 

Table ^ displays results only for when m' = m = 0. To illustrate some of the dependence 
of SCi£rmm' ou the iudiccs m and m', in Figure ^ we plot all of the non-zero entries of 5Cu'mm' 
for £' = £ = ?, and £' = £ = 6. The qualitative manifestation of statistical anisotropy is evident: 
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at fixed t' = £, the correction 6Cwmm' increases with increasing magnitude of m' = m. The 
effect is on the same order as the m' = m = correction to ^cf^)^^,. Recah that in standard, 
isotropic infiation the multipole correlator is independent of m and m'. 

Note that because the elements of dCuimm' do not fall off very rapidly with increasing i, 
there is hope to gather sufficient statistics to detect them, despite their being suppressed by 
roughly the size of the present-day curvature parameter; see e.g. |36]. While we consider it 
very interesting to explore what are the precise limits placed by cosmic variance on detecting 
these effects, such an investigation is beyond the scope of the present work. 



5. Plausibility 

The scenario described in this paper involves the convergence of a number of hypothetical 
circumstances. Yet many of these ingredients are not crucial to the general idea, which we 
expect to cover a broad set of models. Furthermore, while the basic framework — bubble 
nucleation via modulus destabilization followed by not much more inflation than is neces- 
sary to flatten our Hubble volume — is certainly speculative, we do not consider it extremely 
implausible. We here take a moment to justify these attitudes. 

An essential ingredient of our model is that the vacuum in which our bubble nucleates has 
one fewer large spatial dimensions than that of our universe. A toy model to implement this 



is described in Section 2.1; however we consider this model only as a simple illustration of the 
concept. If we accept ten- or eleven-dimensional string theory as a fundamental description 
of nature, then the existence of our (3-|-l)-dimensional pocket is proof of the compactification 
principle, and it seems reasonable to presume that other compactifications are possible and 
that quantum transitions among them occur during eternal inflation. 

The mere occurrence of such transitions does not imply that one is likely to be the 
progenitor of our universe. One might argue that, all else being equal, the nucleation of our 
bubble in a vacuum with a greater number of large (expanding) spatial dimensions is more 
likely, because of the greater physical volume available to nucleate a bubble. However, this 
statement involves an assumption about the spacetime measure on the multiverse, raising an 
unresolved issue in the understanding of eternal inflation. We now briefly explain. 

One consequence of eternal inflation is that the number of bubble nucleations (or of any 
type of event) diverges with time. Indeed, the physical three-volume on an FRW foliation 
within a given bubble diverges as well. Attempts to regulate these divergences have revealed 
that cosmological predictions tend to depend on the choice of regulator — this is known as 
the measure problem (for some reviews see e.g. |]57| , ^). While it is unclear what is the 
correct measure on the multiverse, certain measures can be ruled out for making predictions 
in wild disagreement with our observations. In particular, measures that grant greater weight 
according to larger inflationary expansion factors suffer from the "Q-catastrophe" (and "G- 
catastrophe" ) (5^, |6^, predicting the amplitude of primordial density perturbations and 
the gravitational constant to be in stark conflict with what is observed. Thus, one should be 
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skeptical of arguments about the frequency of events in the multiverse based only on naive 
comparisons of inflationary expansion factors. 

In fact, the question of how the spacetime measure should address transdimensional tun- 
neling in the multiverse has not yet been addressed in the literature. To support the attitude 
expressed above, we here briefly speculate about three of the leading measure proposals. One 



of these is the comoving probability measure |62, p3], which can be seen as weighting events 



according to the frequency at which they are encountered by the future histories of a given 
worldline. A simple generalization is to specify the worldline including its position in any 
initial compact dimensions, in which case this measure would not grant any additional weight 
to parent vacua simply because they have a greater number of expanding dimensions. The 
generalization of the causal patch measure [64| is less clear, but since it also counts events 
according to their proximity to the future histories of a given worldline (specifically, whether 
or not they reside in the surrounding causal patch), one might guess that it would not give 
a very different result. The scale-factor cutoff measure |^, |6^ counts events only if they 
reside in the finite volume between an initial spacelike hypersurface and a later hypersurface 
determined by a fixed amount of expansion. One way to extend this measure to the case 
of transdimensional tunneling is to track the density of a fiducial "dust" of test particles 
scattered over the initial hypersurface, including over any compact dimensions, defining the 
later hypersurface according to when the density drops below a pre-specified value. In this 
case, the extra volume due to an additional expanding dimension would be canceled by the 
extra dilution of the dust due to that expansion, and again parent vacua receive no additional 
weight simply due to having a greater number of expanding spatial dimensions. 

It is not enough to have an effectively (2-|-l)-dimensional parent vacuum nucleate a bubble 
of our vacuum phase — we also require an appropriate period of slow-roll inflation subsequent 
to bubble nucleation. The analysis of this paper introduces a number of assumptions about 
this period of inflation, but we consider all but one to be simple matters of convenience. For 
example, we assumed that the inflaton energy density is precisely constant during inflation, 
and that density perturbations are sourced by fluctuations in a second, subdominant field. 
Yet we expect these choices to affect only the amplitude and tilt of the resulting spectrum of 
perturbations, both of which are model-dependent parameters that can presumably be tuned 
to match observation by simply picking an appropriate implementation of inflation. 

We also assumed that the inflaton energy density dominates over the modulus effective 
potential immediately after bubble nucleation. This was essential to obtaining a simple 
analytic solution for the background metric during inflation, but with hindsight we can see 
this too is not an important assumption. The balance of contributions to the energy density 
in the bubble affects only the time dependence of the scale factors o and b, which in turn 
affects only the power spectrum Pqg (not the anisotropic "Fourier" mode functions Uqrs)- The 
effect of anisotropy in the power spectrum was found to be subdominant to that coming from 
the mode functions; indeed the anisotropy of the power spectrum is ignored in the results 
of Section 4^, since this allowed us to introduce an additional approximation to speed the 
numerical integration given our limited computational resources. 
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The crucial assumption is that slow-roh inflation lasts long enough to conform to present 
observation, but not so long as to push all of the effects of initial anisotropy outside of our 
present horizon. This requires a sort of "fine-tuning" between the inflaton potential and the 
present age of the universe, which, at first glance, seems to involve an unusual coincidence. 
However, such a "coincidence" may find an explanation in landscape cosmology. 

If indeed our pocket is one among a diverging set in an eternally expanding multiverse, 
then we must be careful to account for all of the selection effects that modulate the likelihood 
for us to make a given observation. These selection effects are ultimately encapsulated in 
the choice of spacetime measure; however in the present case it turns out that all three 
measures mentioned above give similar predictions |6^, Roughly speaking, there 

is a factor — called the "prior" — which gives the relative probability that a random spacetime 
point resides in a bubble that undergoes N e-folds of slow-roll inflation. This is convoluted 
with a second factor — called the "anthropic factor" — which is proportional to the density of 
(appropriately defined) observers in bubbles characterized by e-folds of inflation. 

The prior distribution is addressed in ||6^, where an (admittedly crude) argument is given 
to suggest that the landscape might strongly prefer fewer e-folds of slow-roll inflation, with the 
distribution dP{N) tx N~'^ dN. The effect of anthropic selection is most carefully computed 
in [^9| , where it is found that the mass function of Milky- Way mass galaxies is approximately 
independent of l^curv for l^curv ^ C(O.l), and falls off exponentially for roughly Ocurv > 0.7 
(note that Ocurv oc e~^^ with a proportionality constant that depends on the scale of inflation 
and some details of reheating and the subsequent big bang evolution). The combination of 
these effects places most observers in spacetime regions in which ^cuw unobservably small 
(i.e. below cosmic variance), but at the same time gives a reasonable probability, roughly of 
order ten percent chance, to observe ^cuw ^ 10~^ [38, 69]. 



6. Discussion 



The present understanding of string theory and inflationary cosmology points to a picture 
of spacetime containing countless bubbles endlessly nucleating within an eternally-inflating 
background. In the context of the string landscape, the complete set of bubbles contains a 
wide range of low-energy physics, including different numbers of compact dimensions (in fact it 
is the enormous variety of compactifications that inspires the diversity of the landscape in the 
first place). It is therefore a priori possible that our bubble, containing three large (expanding) 
spatial dimensions, nucleated within a vacuum containing only two such dimensions. 

For instance, the eternally-inflating vacuum in which our bubble nucleates could contain 
a compact dimension, the size of which is governed by a metastable modulus that tunnels 
through a potential barrier upon bubble nucleation. The tunneling instanton and initial 
bubble geometry then respect reduced symmetry from the (3-|-l)-dimensional perspective, 
due to the presence of the additional compact dimension. Yet, as expected, a round of 
slow-roll inflation within the bubble is sufficient to redshift away the initial anisotropy and 
curvature, creating an 0(3)-symmetric FRW cosmology consistent with the observed universe. 
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Nevertheless, if inflation within the bubble does not last too long, effects of the initial 
anisotropy may be observable. We here focus on one such effect: the generation of statistical 
anisotropy among large-scale inflationary perturbations. We find that, when projected onto a 
two-sphere approximating the surface of last scattering, the inflationary spectrum generates a 
multipole correlator Cee'mm' that features (in an appropriate coordinate system) off-diagonal 
elements in i and i' (when i — £' = ±2), as well as dependence on the multipole moment m 
(it is still diagonal in m and m'). These effects are suppressed relative to the statistically- 
isotropic components of Cwmm.' by the present-day curvature parameter r^curv; but appear 
to extend to arbitrarily large i. 

There are a number of remaining issues to be explored. Most importantly, as a first 
approach to the problem we have ignored the effects of spatial curvature and expansion 
anisotropy on the free streaming of photons from the surface of last scattering to the point of 
present detection. In fact, anisotropic spatial curvature sources anisotropic expansion, which 
in turn deforms the surface of last scattering away from the surface on which we project the 
inflationary spectrum, in addition to perturbing the trajectories of geodesies as they radiate 
away from the point of observation. A full understanding of the observable signatures of 
anisotropic bubble nucleation requires combining both of these effects. 

Also, we have ignored metric perturbations, focusing on the spectrum of a subdominant 
scalar field and assuming its isocurvature perturbations translate directly into adiabatic den- 
sity perturbations. While we do not expect this to have a large effect on the spectrum of 
statistical anisotropics (compared to for instance the standard scenario where the primordial 
perturbations are sourced by the inflaton itself), it does not allow us to study the tensor 
perturbations generated during inflation. Because scalar and tensor metric perturbations in 
general do not decouple in an anisotropic background, there are possibly interesting correla- 
tions between these signals. Ignoring metric perturbations also hides any interesting effects 
that might come from fluctuations in the bubble wall itself. 

It would also be interesting to explore the nature and degree of non-Gaussianity implied 
by the existence of statistical anisotropy among inflationary perturbations. 

Finally, another potential signature of multiverse cosmology is observable bubble colli- 
sions, see e.g. 0, |72|, [T^]. It would be interesting to understand whether the reduced 
symmetry at early times of anisotropic bubble nucleation affects the spectrum of bubble col- 
lisions on typical observer's sky, or if there is any special signature of collisions with bubbles 
containing a reduced number of large spatial dimensions |76]. 
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Note added: We very recently became aware of interesting work by another group, which 
computes the effect of late-time spatial curvature and expansion anisotropy on photon free- 
streaming from the surface of last scattering [ 77 1 . 



A. 4d anisotropic tunneling instanton 

We here verify the description of the tunneling instanton in the 3d effective theory of Section 
|2.2| , by solving for the corresponding solution in the full, anisotropic 4d geometry. Our starting 
point is the 4d action, (|2.1D. The 4d line element of the parent vacuum can be written 



ds^ = dcj2 + a^{a)[-dt^ + cosh2(t) d(/>2] + b'^{a) dz"^ , (A.l) 

and as before we study the approximate complex-scalar- field solution ip = r/e*"^. Taking 
account of the energy momentum tensor of this model, 

T,, = K'd,^*d,^ + g,, [-Ik{X) - ^ - r?2)2 - , (A.2) 

where X = d^<f*d^ip = •n?"rf'b'''^ and K' = dK/dX, we obtain the equations of motion, 

Above we have used K{X) = X + K2^^ + f^?,^^ as in the main text. These equations of 
motion permit a solution of the form 

a{a) = Hp'' sm{Hp a) , b{a) = bp , (A.5) 

with Hp and bp being constants, provided that we impose the condition, 



S-kG 52 266 
the solution of which is 6^;. The constant Hp is then given by 

ff? = 8-G(^ + ^ + ^)+A. (A.7) 

In the language of the 3d effective theory, this condition can be written v' = therefore 
we have identified the compactification solution but from the vantage of the full 4d theory. 
As expected, inserting this solution into ( A.l ) gives the line element of (2-1-1 )-dimensional de 
Sitter space crossed with a fixed-circumference circle. 
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Figure 4: The 4d instanton solution for a{a) (left panel) and 5(cr) (right panel), all quantities given 
in units of G. The large numbers are due to the dynamics being sub-Planckian; see text for details. 



Consider evolving the equations of motion numerically, taking the boundary conditions 

„(<„ = „ + 2.0(^ + 4^ + ^-— + (A,8) 

,M = ,,_2,g(-^ + ^ + ^-:^).^ + ,,,. (A.9) 

in the limit a — )■ 0. This corresponds to a Taylor expansion of 0(0") and 6(ct), taking 6(0) = 6^, 
and obtaining the other coefficients of the expansion by inserting into the equations of motion. 
Given these boundary conditions, one can show the geometry is smooth in the limit cr — t- 0. 
We then set the value of by trial and error, so that the entire solution is smooth, in 
particular so that as a approaches some value cJmaxi a — ^ Cmax — cr and h — )• constant. The 
results of such a numerical evolution are displayed in Figure ^. To generate these curves, we 
have used the same values of parameters as are used to generate Figures Q and ^ in Section 



The results agree with those of Section in the sense that when we read off the initial and 
final values of the circumference of the z dimension, between the two approaches they agree. 

While the interpretation of this instanton as a tunneling event is not so apparent in the 
4d picture, note that the same instanton would describe the reverse "tunneling" process, 
the description of which is less clear in the 3d effective theory, which breaks down inside 
the nucleating bubble. (While we do not include the evolution here, there is no difficulty in 
extending the evolution of a and h to convey the evolution in the daughter vacuum.) Viewing 
this solution as a tunneling instanton, one can interpret it as an inflating black brane, charged 
with respect to the scalar field ip, nucleating in the 4d de Sitter background spacetime. Higher 



dimensional solutions analogous to this have been discussed by [ |12| . 
B. Bunch— Davies vacuum 

To determine the power spectrum of the light scalar a with respect to the Bunch-Davies 



vacuum |52], we trace the evolution of T^^ back into the parent vacuum and Klein-Gordon 



normalize the positive- frequency modes over a Cauchy surface, as is done in |52, 54, 55]. 
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Figure 5: Conformal diagram of a de Sitter bubble nucleating in a de Sitter parent vacuum, indicating 
the bubble wall (thick solid curve), the future lightcone of the bubble nucleation event (thick dashed 
curve), some surfaces of constant rjc and rj^ (solid curves), and some surfaces of constant x x 
(dashed curves). The uncharted regions of spacetime are not relevant to our discussion. The Cauchy 
surface x = is indicated. 



Figure H displays a conformal diagram illustrating the geometry (the coordinates z and (/) 
are suppressed). For simplicity we have constructed the diagram as if the vacuum energies 
inside and outside of the bubble are the same; however the statements below apply equally 
to more realistic geometries. We have also drawn the diagram (and we work below) as if all 
of spacetime can be covered by just these two universes — in other words, as if there were no 
other bubbles.^ Within the bubble, which nucleates at the middle of the left boundary and 
thus occupies the upper left corner of the diagram, x is a radial coordinate on constant-z 
(2+l)-dimensional open FRW hyper surf aces, and ijc provides a spacelike foliation. Evidently 
no Cauchy surface can be drawn within the bubble. 

To simplify the geometry of the parent vacuum we assume the energy density of the 
bubble wall does not significantly perturb the metric. Furthermore, we treat the radius bd 
of the compact dimension z as if it were essentially constant throughout the parent vacuum, 
and we take the energy density immediately inside the bubble to be essentially equal to that 
of the parent vacuum on the outside (to be precise, we take H?^^ = A/3 to be essentially 
the same as the Hp of ( |A.7D ). Note that these are not necessarily good approximations of 
the toy compactification model described in Section ^; however they greatly simplify the 
computations here. In particular they allow us to treat the parent vacuum as pure de Sitter 
space (crossed with a circle), in which case the central, diamond-shaped region of Figure ^ 

*In fact the daughter bubble experiences an infinite number of collisions with other bubbles that nucleate in 
the parent vacuum ^ |7^, |7^ ; however the phenomenological success of inflationary theory suggests that 
any effects of such collisions should be small — linear perturbations to the standard inflationary background (at 
least in our past lightcone). Meanwhile, the parent vacuum itself likely resides in an open bubble, or contains 
other bubbles that intersect the Cauchy surface at x = (see text below). However bubble nucleation rates are 
likely to be exponentially suppressed, so any such disturbances would typically be very far from the physical 
region of interest. For simplicity we here presume both effects are negligible. 
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can be covered by coordinates with line element 

ds'^ = a2(r/J [diql - dx'^ + cosh^{x) +bldz^. (B.l) 

Here the coordinate x provides the spacelike foliation, rj^ is a spatial coordinate, and the 
hypersurface x = is a Cauchy surface. Since x = lies entirely within the region covered by 
( B.l ), it is unnecessary to cover any more of the spacetime by introducing additional charts. 

Note that the bubble coordinate system of (O) can be smoothly connected to the parent 
vacuum coordinate system of (B.l) by using analytic continuation, 

— - iir _ . ,^ 

X = X-y, Vc = Vc + — , a = ia. (B.2) 

Here the use of analytic continuation is merely a technical device to describe the propagation 
of modes through the bubble wall. In particular, we do not transform the coordinate of the 
bubble scale factor b{r]c). Instead the z dimension of the parent vacuum maps directly onto 
that of the bubble (as is the case with the coordinate 0), assuming that b{r]c) — >• bd and 6^ — )■ 
as r]c —5- — oo. 

Within the parent vacuum we denote the mode expansion of a according to 

'a bd V« bd 

where to simplify the presentation here and below we do not worry about keeping track of the 
(unimportant) overall phase. We have introduced the superfluous normalization factor Mq 
for later convenience. The (p- and z-dependent Fourier modes are the same as those within 
the bubble, $r(</') = ^r{(p) and Zs{z) = Zs{z). Meanwhile, the x-dependent modes can still 
be written as solutions to the Legendre equation, but they now take the form 

^^^^^^ = ^^"^-g)"'^ ^i^^(^)) • (B.4) 

The normalizations of ^r-, Zg, and Xqr are such that, after analytic continuation, each mode 
maps onto the similarly-denoted mode defined within the bubble. Note that the modes Xqr, 
which describe the time-like evolution in the parent vacuum, correspond to positive-frequency 
plane waves oc e"*''-^ (up to an unimportant phase) in the limit x ~^ 0. 

The analytic continuation of the scale factor a gives a = sech(r/(,), while in the parent 
vacuum the scale factor b is constant, b = bd- This gives the equation of motion 



CTvqrsiVc, X, 0, z) = -^== T^^ (r? J C/gr-s (x) = — ^ T^s(^c) ^gr-(x) ^r(</>) Z s{z) , (B.3) 



T + 

^ qs ~ 



T:, = 0, (B.5) 



where the dots now denote derivatives with respect to fJc, and again = s/bdH\^i. Notice 
that the separation constant /x'^ appears as would a mass for the scalar field a. The general 
solution of (B^) can be written in terms of associated Legendre functions: 



C\Fi{r),) + C"^F2{ri,). (B.7) 
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The solutions Ti and T2 represent two orthogonal spatial modes, and thus each should 
be normalized separately, with the two summed over in the final mode expansion. Hence we 
have introduced the index u, which distinguishes these modes. To perform the normalization 
it is convenient to apply a trick. Notice that ( |B.5| ) has the form of a Schrodinger equation for 
a point particle with position r]^, potential energy = - 3/4)sech2(7yJ, and energy 

q^. The potential Virj^) tends to zero as r/^ — >■ ±cxd; thus T^^ may just as well be expressed 
in terms of two linearly independent "scattering" solutions of the form 



iT 



qs 



iT 



e 



iqVc 



qs 



UJ+ e*"?^- + e- 



as r/^ 

as fj^ 

as r/^ 

as r/^ 
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+00 , 

— 00 , 
+00 , 



(B.c 



where studying the Wronskian reveals |w±P + |?±P = 1, = <^_, and <j+a;l = |55]. 

(The factor J\fg was introduced in (4.4) and (B.3) to permit these simple normalizations; it 
will be determined later.) The normalization of the delta function is determined by noting 
that all of its support comes from large \f]^\, where the solutions are plane-wave. Thus 



dVc'^L%*s = 27r 6{q - q') 5„ 



(B.IO) 



where the index v takes values it. We have introduced the modes T^^ merely to take advantage 
of this simple normalization. 



The normalized scattering solutions T^^ can be expressed in terms of the solutions Ti 



and J-2 by matching their asymptotic behavior onto (B.8) and 
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By matching the asymptotic positive and negative "frequency" modes of Ti and J- 2 onto the 
asymptotic behavior of T^^ and T^^ given in ( [B.8D and (p3.9| ), we determine 



c4 P2 - "2" /^r 



T,-.(5j,) = (B.17) 



(Some of the terms above are zero, however we delay simphfication until later.) 

We have now determined all of the mode functions, up to the overall normalization Mq. 
This is determined by enforcing Klein-Gordon normalization of the scalar field o". 



' v'q'r' s' \ 



iPvqrsi O'v'q'r's') ~ ~^ j ^'^P'-d^'^ \_<^vqTS du<^v' q'r' s' ~ i^uO'vqrs) C^v' 

= ^{Q - q') ^rr' Sss' 5vv' ■ (B.19) 

As mentioned before, we set the Cauchy hypersurface at x = 0. Using the orthogonality 
of the mode functions T^^, and Zg, the above normalization condition can be written 

27r \Mq\^ COsh(x) \Xqr B^^r " i^'^ir) ^r] _ = ^ " (^-20) 

The term in brackets is at first glance rather complicated; however it can be simplified with 
some technical manipulations, and in the end we find (up to an unimportant phase) 

Mq = ^ . (B.21) 

Y'4g sinh(7rg) 

Now that we have the Bunch-Davies modes of a in the parent vacuum, the next step 
is to propagate these modes into the bubble. For the modes and Zg this is easy: they 
are unchanged. The modes Xqr and T^^ are propagated by analytic continuation of the 
coordinates in their arguments, as given by (|B.2D . The temporal mode functions Xqr{x) 
become the bubble spatial modes Xqr{x) of (|4.8|) . It is left to discuss the modes T^^. 

Note that the only dependence of the T^^ on r]^ is via tanh(?7g), which is rotated into 
coth(7/c) when we take 'nc + ivr/2. The functions J-'i and 7^2 thus become 

^i(f?,) ^ ^1(7?,) = P^^_^^^(coth(r?,)) (B.22) 
T2ir]c) ^ Mric) = ^ , (coth(r?e)) . (B.23) 

These are not the same as the functions J-i and J-2 computed with the bubble line element 
([4. 3D, because those solutions account for the growth in the scale factor b{rjc), whereas Ti and 
J-2 continue from the parent vacuum where b is static. Nevertheless they have the correct 
asymptotic form, as tJc — )• —00, because the scale factor 6(r?c) oc coth(r/c) approaches a constant 
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in that limit. All we must do is match the solutions T\ and T2 onto T\ and J-2, in the limit 
r^c — >• —00. In this limit, T\ and T2 have the asymptotic behavior 
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+ /3i e 



where the coefficients dj and /3j are given by 
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(B.25) 
(B.26) 
(B.27) 
(B.28) 



The functions J-\ and J-2 have been expressed so as to have simple asymptotic behavior, 
see ( [1.17 ). Putting all of this together, we obtain the pair of mode functions 
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(B.30) 



Referring to (p35|)-( p38|) , we see that the terms in the first parantheses of ( |B.29| ) cancel to 
zero, as do the terms in the second parantheses of ( [B.30| ). Referring also to ( [B.12| )- ([B.16 ), 
we see that other basic simplifications are possible, and in the end we obtain 



(B.31) 
(B.32) 



The power spectrum can now be computed just as in Section |4.1| , except now we must 
include the extra normalization factor J\fq. That is, we write 



{a{r]c, q, r, s) a\r]c, q' , r', s')) 



■ Yl - 1') ^rr' 6ss' ■ (B.33) 



As described in Section 4.1 we evaluate the power spectrum in the limit r]c — >• 0, in which case 
the asymptotic behavior of the time-dependent parts is given by and (|4T9| ). Putting 

everything together, and performing some manipulations to eliminate some of the gamma 
functions, we find the power spectrum to be 
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Figure 6: The power spectrum P^s, expressed as hY{hdPqsl Hi^{), plotted as a function of g for /i = 0.1 
(top pair of curves), /i = 2 (middle pair of curves), and /i = 10 (bottom pair of curves). In each case 
we plot the actual power spectrum (solid) and the "asymptotic limit" referred to in the text (dashed). 
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In Figure ^ we plot the above power spectrum, and its small scale (large q, ^) asymptotic 
limit ( [4.26 ), for a few values of ji. It is clear that the spectrum rapidly approaches its asymp- 
totic limit, with significant deviations only for q < 0(1). This is not unlike the situation in 
regular (isotropic) 4d open bubble inflation, where the inflationary spectrum is also modified 
s-t 9 ^ C'(l) ||5^, |5^, 55]. In the present situation, it is evident that the deviation from the 
asymptotic curve is itself anisotropic; for example at g = 2 the curve with ^ = 0.1 deviates 
more significantly than the curve with = 10. This anisotropy is on top of that which results 
from the use of anisotropic mode functions in the power spectrum. 
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